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Starting from classic work of Feynman on the A-point of liquid Helium, we show that his idea of 
universal action per particle at the BEC transition point is much more robust that it was known 
before. Using a simple "moving string model" for supercurrent and calculating the action, both 
semiclassically and numerically, we show that the critical action is the same for noninteracting and 
strongly interacting systems such as liquid ''He. Inversely, one can obtain accurate dependence of 
critical temperature on density: one important consequence is that high density (solid) He cannot 
be a BEC state of He atoms, with upper density accurately matching the observations. We then use 
this model for the deconfinement phase transition of QCD-like gauge theories, treated as BEC of 
(color)magnetic monopoles. We start with Feynman-like approach without interaction, estimating 
the monopole mass at Tc. Then we include monopole's Coulomb repulsion, and formulate a relation 
between the mass, density and coupling which should be fulfilled at the deconfinement point. We end 
up proposing various ways to test on the lattice whether it is indeed the BEC point for monopoles. 

PACS numbers: 



I. INTRODUCTION 

The goals of this paper are two-fold. The first goal 
is rather general: to get better qualitative understand- 
ing of the parameters controlling the transition between 
the "normal" matter and its Bose condensed versions, 
for strongly interacting bosons which may be in form of 
a liquid (fluid-superfluid transition well known for *He) 
or solid (solid-supersolid transition yet to be found) . Al- 
though there are high quality Monte Carlo numerical re- 
sults for ^He and many other systems, we think the uni- 
versal condensation criterion is still very much needed, 
as performing numerical simulations is not trivial in each 
new settings. 

As it will be explained in detail below, we will fol- 
low 50-year-old Feynman theory of Bose condensation 
[U [2] , in which he introduced the notion of the critical 
value for the jump amplitude or the critical action 
Uc = exp{—Sc)- When Feynman realized that his sim- 
ple treatment (evaluation of only the kinetic energy part 
of the action) needs correction, he simply introduced an 
"effective mass" thinking that some extra matter is in- 
corporated into exchange motion. We think instead, that 
other particles (except the ones in the exchanged poly- 
gon) have very little chance to move. Instead, the jump 
amplitude should be correctly evaluated, with the inter- 
action term included. We thus revive Feynman's idea, 
using a simple model of particle motion -the "moving 
string model" - which can be studied either semiclassi- 
cally or numerically of the smallest-action paths which 
particles should follow during their exchanges. 

Our second goal is very far from atomic systems: it is 
related with the deconfinement phase transition in QCD 
and related gauge theories. "Dual superconductivity" 
or BEC of certain magnetic objects in the vacuum of 



these theories were proposed to be responsible for con- 
finement by t'Hooft and Mandelstam [H]. As shown in 
refs f5^, TT] , confining strings (electric flux tubes) can be 
well described by effective models, making them dual ver- 
sions of the Abrikosov flux tubes in superconductors [6]. 



As Feynman did in 1950's, we would however approach 
the problem from the "normal" phase above the decon- 
flnemcnt phase transition called "Quark Gluon Plasma" 
(QGP). At very high T one can view QGP as a plasma 
made of "electric" quasiparticles, quarks and gluons. 
However at lower T it becomes a "dual plasma" contain- 
ing not only electrically charged quasiparticles but also 
magnetically charged objects - monopoles and dyons. 
As it was pointed out in [TU], with a decrease of the 
temperature and increase of the electric coupling con- 
stant one expects a gradual shift from electric to mag- 
netic dominance, with large density of (color-magnetic) 
monopoles near Tc- Recent lattice studies, in particularly 
[12 , have discovered many important details about such 
monopoles. They indeed found rather high densities of 
them close to Tc, as well as clustering behavior similar to 
BEC. They also measured equal-time density correlators 
of like and unlike monopoles (Fig. 5 of [T^ ) which clearly 
shows the peaks characteristic for strongly coupled ionic 
liquids. The estimated magnetic Coulomb coupling does 
indeed show running opposite to electric coupling pre- 
dicted in ^lOj: see review [7J. So, as QCD deconfinement 
transition is interpreted as BEC of monopoles, in the sec- 
ond part of the paper we will apply Feynman's universal 
action idea to constrain the properties of these monopoles 
(mass, density and coupling constant) for BEC to become 
possible at the deconfinement point. 
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FIG. 1: Example of polygon due to particle exchange 



II. REVIVING THE FEYNMAN'S THEORY OF 
THE HELIUM A-POINT, 55 YEARS LATER 

A. Feynman's theory 

Feynman's starting point is the thermal partition func- 
tion in the form 
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where m* is the effective mass of a He atom, the sum 
is done over permutation P of the particle coordinates, 
the function p includes effects of the interparticle interac- 
tions and the exponent contains additional action due to 
kinetic energy. His main idea is that the function p can be 
inferred from general properties of the liquid, its quasi- 
ordered local structure with peaked distribution over in- 
terparticle distances at some nearest-neighbor value d (in 
the case of cubic lattice d is the lattice spacing). Then 
the relative magnitude of a term with permutation of n 
atoms would be proportional to the n-th. power of the 
"jump amplitude" which is approximated as 



Uf = exp 



ra*Td'^ 
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with some combinatorial prefactors, describing a num- 
ber of corresponding to non-crossing polygons in the 3d 
lattice. The divergence of this sum, at the parameter 
approaching some critical value uf ^ Vc should indicate 
the presence of "infinite cluster", the signature of Bose 
condensation. The condition for this, Feynman argues, 
is He ~ 1/s where s is the growth factor in the number 
of polygons when n is increased by one unit. Feynman 
mentioned that he expected ?/c = l/4-l/3. 

The combinatoric problem to find the correct critical 
value of the Feynman parameter, was studied in details 
by Kikuchi et al. [J. They found a critical action Sc ~ 
1.9. We can consider the distance d of a "jump" fixed 
to the position of the nearest neighbor maximum in the 
static correlation function g(r) for liquid He, which is 



d « 3.5 A. If this value is used and the critical action 
is taken to be the one obtained by Kikuchi, in order to 
recover the correct position of the A-point Tc = 2.17, 
the Helium mass is changed to an effective value m* — 
1.64mHc- Otherwise, using the physical Helium mass we 
obtain Tc — 3.57 which is, of course 64% bigger than 
the real critical temperature. In the eighties Elser [3] 
determined numerically that the critical action should be 
smaller than the values obtained in S";, « 1.44. Using 
this new parameter the predicted critical temperature 
reduces to = 2.72, closest to the physical temperature. 

In the Feynman picture all the effects of the potential 
are absorbed into an effective mass m* that appears in 
the kinetic part of the partition function. In this letter 
we try to study the He A-point, starting from the idea of 
Feynman, but including the effects due to the potential 
at a mean field level. In our calculations we have con- 
sidered the Aziz potential HFDHE2 [5]. We will show 
which in this context an estimate for the Helium criti- 
cal temperature can be obtained within an error of 5% if 
we assume that the critical action will be the same than 
the free case, that, as we will show, can be estimated to 
Sc = 1.655 for a cubic lattice. 

Following Feynman we put the Helium atoms filling a 
cubic lattice and consider two different states for the sys- 
tem: the first where the atoms are static on the lattice 
sites, the second where the atoms on a line of the lattice 
are moving coherently in a given direction, for example 
along the x, see Fig. [2] If we consider an infinite sys- 
tem we can consider that this line is a side of a polygon 
with infinite number of atoms. Another difference with 
the Feynman's work is that in our case the atoms lying 
on the moving line are not restricted to move on the lat- 
tice link: what they must do is to jump from one site 
to the following on the same line, but in between they 
move following the three-dimensional equation of motion 
imposed by the geometry of the selected configuration. 

The appearance of such a kind of particle cluster, as 
proposed by Feynman, is connected to the transition into 
the superfiuid phase. Confirmation of this hypothesis 
come also from Path Integral Monte Carlo (PIMC) cal- 
culations, which, at today, represent the only way to de- 
termine without any approximation, the Helium critical 
temperature. In this framework, the transition point is 
fixed exactly from the appearance of atoms involved in 
an exchange which winds around the simulation box. 

The idea of Feynman was to consider the ratio between 
the partition function of the system where permutations 
are possible with respect to the " Boltzmann" case where 
particles do not permute. In our case we are interested to 
understand the properties of the particular system where 
there is exactly one line involved in permutation with 
respect to the Boltzmann system. 

Systems at finite temperature can be studied using 
path integral techniques in the Euclidean space, where 

t —IT. 

In general the partition function for bosonic interacting 
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FIG. 2: 2 = of the three dimensional lattice: the red con- 
tinuous arrows identify the moving atoms. The dashed arrows 
represents the diagonal trajectory of monopoles discussed in 
Sec.lml 

particles in Euclidean space can be written as 
...xn--,P) = 

P 

where [3 — and the Euclidean action is 

SE[x„..,x^;r]^J2{^-i(-) + T.^^^)^ (4) 

i=i j=i 

and only two-body interactions are taken into account. 

Now we consider the case where K particles 
(xfc, ...Xfe+x) on a line move coherently, that means: 
Xk+jir) = Xfc(r) -I- jd. This hypothesis implies 
that we have chosen the particular permutation where 

Xi^k,....k+K{T) = Xi,Xk{T) = Xk+l{0), .■.,Xk{t) = Xfe(O). 

Moreover we impose that the position of all other parti- 
cles is constant. 



B. The "moving string" model 

Having selected this particular configuration, we can 
evaluate the potential acting on a particle lying on the 
line, due to all the other particles. 

In Fig. [3] we can see the projection on the x — y plane of 
the resulting mean field potential and one of the possible 
path followed by the atoms on the moving line. 

In particular we consider only the a;-coordinate and 
study the dependence from the density of the mean field 
potential. In Fig. [4] we can see how change the behavior 
of the potential considering the atoms disposed on a cubic 
lattice and changing the density. In Fig. [5] we do the same 
for an HCP crystal. We can see that there is a range 
of density in both case where the origin is no longer a 
minimum, which means that the lattice considered is not 
the correct one to study the system. Therefore we need to 
pay attention which configuration we choose depending 
on the density. 




FIG. 3: Projection on the x — y plane of the potential for a 
particle on the moving line 
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FIG. 4: Density dependence of the potential when we con- 
sider the atoms on a cubic lattice 
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FIG. 5: Density dependence of the potential when we con- 
sider the atoms on a HCP lattice 



Another fact related to the Aziz potential is repre- 
sented in Fig. [6] where we see that the amplitude of the 
sinusoidal potential does not increase continuously when 
the density grows. We can think that such a behavior 
will be connected wit a wrong description of the atoms 
distribution in this range of densities. 
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FIG. 6: Density dependence of the amplitude of the sinu- 
soidal potential. 
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FIG. 7: Density dependence of the NP potential for five 
different values of p. The point can be successfully fitted by 
the function Vo(n) = anF^'^ . 



This is indeed peculiar of the Aziz potential where we 
have both attractive and repulsive terms. If for example 
we consider a simple negative power potential, such a 
behavior not only disappears but we can find a simple 
relation between the amplitude Vq and the density. 

Define the negative power (NP) potential 
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with Op > 0. The scattering length can be defined for 
this potential when p > 3 and is given by 



/2map/^\ '/'^-'^ r[(p-3)/(p-2)] 
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where T[x\ is the Gamma function. We consider p = 
1, 4, 6, 9, 12 and choose to keep fixed the value a = 1 for 
all the potentials, changing accordingly ap. For p = 1 we 
arbitrarily fixed ai = a^. For this potentials we compute 
Vb for different density, taking care that for such a NP 
potential the atom disposition must be a bcc crystal. The 
result are presented in Fig. |7] the data can be fitted with 
extraordinary accuracy using the function: 



an 



p/3 



(7) 



We can try to use the same function also to fit the 
sinusoidal amplitude obtained by the Aziz potential, in 
this case we can use Vb(7i) — amP^^ where p = 6 is the 
dominating term for large distances that appear in the 
Aziz. The result is in Fig. |8] we can see that the fit is 
quite good until, as expected, we enter in the yellow band 
described before. 

Nevertheless we decided to use this simple scaling re- 
lation for our calculations. 

Due to the particular configuration we have chosen, the 
partition function for the Bose system can be rewritten 
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FIG. 8: Density dependence of the Aziz potential fitted by 
the function Vo(n) — an^l'^ . 
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where Co{x) is the potential acting on the particle at rest 
and which is exactly the potential that also appears in 
the Boltzmann partition function. 

For this reason the ratio between the Bose (Z) and the 
Boltzmann (Zb) partition function reduces to 
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and the Feynman parameter is given by 
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C. The semiclassical calculation 

The simplest case is that of noninteracting particles, 
for which the optimal path ("Feynman's instanton") is 
just the straight line starting at the initial position of a 
particle and ending at the previous position of its neigh- 
bor (as indicated by lines with arrows in Fig. [2]). Since ve- 
locity is constant, the extra action per particle needed for 
its jump is nothing but just the kinetic energy times the 
Matsubara time (3 = h/T available for the interchange, 
with an obvious velocity on the segment being d/P 
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per particle, for each of the diagrams. Let us see what is 
this action at the Bose-Einstein condensation (BEC)line, 
for the ideal gas. Using 



2/3 



/to 



(12) 



and relating the distance of the jump to the density by 
d — Kn~^^^, with some coefficient k = 0{1) one finds 
that the critical action for the "Feynman instanton" is 



'BBC 



/h- 



Kn.655 



(13) 



The precise value depends on the numerical constant k, 
which in turn depends how exactly particles are corre- 
lated in space. It is unclear what should it be for ideal 
gas, but for strongly coupled systems we are mostly in- 
terested in it is uniquely determined by the type of local 
crystal structure developing, e.g. it is exactly k = 1 for 
cubic crystal. We will use below this value, although 
there is some uncertainty here for weakly coupled sys- 
tems. 

When density is too low or T is too high, so that 
S > Sbec, the sum over the polygons is exponentially 
convergent and thus the gas remains in normal phase. 

Now we switch on the potential and impose our atoms 
to be displaced on a cubic lattice. As a first analytic 
approach to the problem we look for the semiclassical 
tunneling path - known as periodic instanton, or caloron. 
It corresponds to solution of classical equation of motion 
in Euclidean time, now including the potential. All cases 
we are interested in this periodic potential on various 
lattices and forces happens to be well described its first 
harmonics, the sinusoidal potential: 



V{x) = VosiT?{xTT/d) 



(14) 



which interpolates between the minima of V 

- 

XciiP) = d 
Introducing conserved Euclidean energy Ee 
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after separation of variable we has the solution 
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where F is the elliptic integral of the first kind. Imposing 
Xci{tq = 0) = we have 
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with JA the Jacobi amplitude for elliptic functions. Par- 
ticular solution for Ee = is called the instanton, it 
corresponds to the zero temperature ( or ^ oo) limit. 
In our setting it is not interesting since we only approach 
the critical point from above and thus T > T^. 

For arbitrary temperature the duration of the Matsub- 
ara time 13 is prescribed: thus from Eq. [17] we obtain the 
equation to fix the energy. One can also read it imme- 
diately as the critical temperature as a function of the 
parameter of our system (d, Vq and Ee) 
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provided the energy is substituted from the Feynman 
condition, which say that the Euclidean action 
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is equal to the Feynman critical value. 

In the instanton limit Ee = we have for the action 
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And the corresponding trajectory is given by 



x{t) = —d arccot 
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In Fig. [9] we can see the difference between the instanton 
trajectory and the simple straight hne corresponding to 
the case of a free Bose gas. 




FIG. 9: Instanton trajectory we have in the presence of a 
potential barrier, compared to the free Bose gas case. 

More detailed comparison between real He and 



caloron/instanton approximation is shown in Fig. 10 'b) 



While it is qualitatively correct, it is not quite accurate. 
This is by no means surprising: it should be accurate 
provided the critical action for BEC would happen to be 
much larger than 1, while unfortunately it is only 1.65 or 
so. Therefore the semiclassical result should be used for 
qualitative comparison only. 

One of those is that as the matter is further compressed 
by extra pressure, its density and the amplitude of the 
potential Vo(n) grow. Respectively the action required 



for tunneling grows, see Fig. 10 'a), and when the action 
gets too large So > Sc, Feynman's condition could not be 
fulfilled. Therefore, a sufficiently compressed "^He cannot 
support a supercurrent of the *He atoms. 

(For clarity: we do not make any statements here 
about possible "supersolid" behavior of the sohd *He in- 
duced by supercurrent of some defects/dislocations im- 
posed on it. All we are saying is that the '^He atoms them- 
selves do not create such supercurrent: this statement is 
however not new and it has been verified in dedicated nu- 
merical Monte-Carlo studies. The only new element in 
our statement is that it follows from Feynman's universal 
action.) 

One may improve the semiclassical expressions for the 
tunneling action in one or two loop quantum corrections. 
For example with the one-loop accuracy the Feynman 
parameter can be written as 



ifi2 2E 



^0 



(detF[xci]) 



1/2 



(23) 



2.0 



1.8 



1.6 



12 



So> 1 .6S5 

cannot condense . 



Sii< 1.655 
can condense 



DIO 0.012 0.014 0.016 0.018 0.020 0.022 0.024 




1 X 10"'' 



5x10-^0.001 



0.005 0.010 



FIG. 10: (a) The action for the tunneling (caloron) solution 
as a function of the density. When the action is larger than 
the critical value S = 1.655 (dashed line) the system cannot 
be Bose-Einstein condensed phase. 

(b) The critical temperature obtain from the 
caloron/instanton solution (dashed line) and from the 
1-d PIMC simulation (points). The red star shows the 
physical location of the lambda point. 



Expressions for two loop are a bit more involved, see 
Ref. [T7]- However we have not calculated it analytically, 
using numerical path integral instead. 



D. The 1-d quantum path integral 

As we have already emphasized above, since the ac- 
tion is not large, one cannot expect the semiclassical the- 
ory be really accurate. Fortunately quantum fluctuation 
around the classical trajectory can be taken into account 
by numerical evaluation of the r.h.s. of ( 10 1 can be sim- 
ply performed using PIMC code. Note that compared 



to PIMC simulation of manybody these calculations are 
very cheap in terms of computational power because we 
consider a one dimensional system with only one particle. 

In order to find the critical temperature as a function 
of the density we do the ansatz that the value of the Feyn- 
man parameter of our system is more or less the same as 
for the free Bose gas. Because we know that for the free 
Bose gas Sp = — log(yF) = 1.655 independent from the 
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FIG. 11: Density dependence of the critical temperature. 
Black points are obtained using the exponential dependence 
of Vb(w).The red star is the physical result for ''He. The blue 
square is obtained using the Vq as obtained directly from the 
Aziz potential. 



density, we simply computed S for different temperature 
at fixed density and when S{T,n) = 1.655 we simply 
called that critical temperature, in Fig. [TT] we plotted 
the results obtained in the case of the cubic lattice. We 
can see the good agreement between our calculation and 
the physical critical point. Also when we use the formula 
Vo(n) = an^/^ for the prediction of the amplitude, our 
prediction is only four percent larger than the physical 
value. 

Is also interesting to notice that the physical point is 
the point where the critical temperature is larger, in this 
sense is looks like natural that condensation happens in 
nature exactly at that point. 

Now we plot the ratio T^/Tq (Fig.[lO](b) black points), 
where Tq is the critical temperature for the free Bose gas. 
In the limit of vanishing density we recover the Einstein's 
critical temperature as it should be. In the region of small 
density there is universal theory using small parameter 
a'^n, with a being the scattering length. As shown e.g. 
in Ref. [TB] , there is a nontrivial increase of Tc by several 
percents there, in agreement with numerical studied. Our 
model is obviously too crude to reproduce it. 

We are mostly interested in the opposite limit of high 
density: here the considered ratio decreases until a cer- 
tain point. For larger density we simply cannot have con- 
densation in the sense that we cannot find any tempera- 
ture where S{T) = 1.655. This can be seen more clearly 



looking at Fig. 
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Starting from around n ~ 0.033 A~'^ 
we cannot find any critical temperature different from 
zero. As a consequence we have a jump from a minimum 
of Tc ~ 1.7 K to Tc = around that density. 

"'He phase diagram (Fig. 



13 1 is also 
at P = 0, 



Looking at the 
interesting to notice that the A-line starts 
T — 2.17 and ends at some pressure with T — 1.76 very 
close to our estimate for the minimal critical tempera- 
ture. Obviously the density where we observe this be- 
havior is completely different, nevertheless is interesting 
to notice that the two minimal critical temperature coin- 



FIG. 12: Action as a function of the temperature for different 
densities. 
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FIG. 13: ■'He phase diagram 



cide. It seems to suggest that it is not possible to have a 
supersolid phase of Helium because, as soon as we enter 
in the non-liquid phase, the critical temperature jumps 
from Tr = 1.76 to zero. 



III. CONFINEMENT AS BOSE-EINSTEIN 
CONDENSATION OF MONOPOLES 

First of all, let us explain why we think Feynman cri- 
terion based on particle paths and clusters is especially 
well suited for this task. One reason is that lattice sim- 
ulations are able to follow monopole's paths, and thus 
identify "single" ones, with individually periodic paths, 
as well as those belonging to k-clusters. Therefore, one 
can see divergence of the cluster expansion at Tc directly, 
without any calculations. 

The second reason is that excited matter we discuss 
has no nonzero quantum numbers, it is just an "excited 
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vacuum" produced in high energy collisions. Thus there 
is no any nonzero conserved charge to which the corre- 
sponding chemical potential can be coupled. Neutral- 
ity leads to equal number of electrically charged quarks 
and antiquarks, as well as equal number of monopoles 
and anti-monopoles. Thus unfortunately one cannot in- 
troduce chemical potential and use the usual reasoning 
related to its crossing the lowest level. 

And yet , since the density and mass of such monopoles 
depend on T, only when the Feynman criterion is satis- 
fied their BEC may happen. In this section we estimate 
what mass and coupling constant the monopole should 
have in order to condense. We do so in two steps, first 
for noninteracting monopoles, including relativistic ac- 
tion, and then for interacting ones using nonrelativistic 
formulae derived above. Again following Feynman, we 
approach the problem from the high-T phase - the QGP 
- in which monopoles are not Bose-condensed, calculate 
the action needed for a "jump" to the site of the identical 
neighbor, and compare it to the universal value Sc- 

We will only consider pure gauge theories (no quarks) 
and ignore gluon quasipaticles. The standard units used 
in this field are based on h ^ c = I, and thus only one 
unit (length in fm or energy in GeV, so fm*GeV=0.1973) 
needs to be defined. It is standard in lattice works to 
calculate the string tension a at zero T and fixed its value 
to be the same as in real QCD, namely ^/a = 0.42 GeV, 
which we will follow. The critical temperature in such 
units is Tc ~ 0.27 GeV. Lattice Tc for pure SU{3) gauge 
theory is in such units about 270 MeV. 



A. Free monopoles 

The work [12] (pure gauge SU(2) theory) has provided 
measurements for monopole (plus antimonopole) den- 
sity n at T = (1 — 12)rc- In this theory Higgsing leaves 
only one massless U(l), and thus there is only one kind of 
monopole. For orientation, near Tc they find n/T^ w 0.3. 
This the distance between identical monopoles is 

a = (n/2)-i/3 = 1 /(0.53 * 0.2706^) « 1.4 fm (24) 

(where 2 is because we only need density of one charge). 

We now estimate the action of the monopole which 
jumps to the position of another identical monopole at 
distance a away, during Euclidean time duration equal to 
the critical Matsubara time (3 = h/Tc = 1/(0.27 GeV) = 
0.73 fm. In the Feynman approximation - when only ki- 
netic part of the action is included - the (Euclidean) ve- 
locity on the optimal path is constant v ~ a/ j3. Putting 
numbers reveals two complications: (i) w > 1 and (ii) rel- 
ativistic monopoles require relativistic form for Euclidean 
action. Since we speak about tunneling, having imag- 
inary action and velocity above that of light is in fact 
appropriate: no negative roots appear. The Euclidean 
action at the BEC Tc point should be 



Putting numbers into the square root one finds 

TO = * S'c/^ ~ 200 MeV (26) 

From existing lattice data on the mass, estimated from 
the paths themselves by D'Elia (private communication) 
one finds rapid decrease of the monopole mass as T ^ 
Tc to similar ball park, but the accuracy of the data is 
not yet sufficient to tell its numerical value at Tc yet. 
And before this comparison is made, we should move on 
into much more involved estimate for strongly interacting 
monopole plasma. 

B. Strongly interacting monopoles 

Color monopoles interact via Coulomb-like magnetic 
forces related to their charges. We remind the reader 
that "Higgsing" is assumed to be due to nonzero back- 
ground field < ^0 >i which lives certain Abelian sub- 
groups of the gauge field unbroken, that is massless. For 
the SU(2) color group there is only one diagonal gener- 
ator which is left massless, so there is only one U{1) 
charge and thus we can use the language of electrodynam- 
ics, simply calling monopoles to be positive and negative 
charges. (This is strictly speaking not true for physi- 
cal SU{'i) color, which has two unbroken Abelian fields, 
proportional to A'^, A'* Gell-Mann matrices, with two fam- 
ilies of monopoles.) One more comment is that in general 
forces include also Higgs exchange, which in principle can 
also be long range if Higgs is massless. However the data 
about monopole correlations [12j seem to show only pure 
Coulomb- like forces. 

Such charges obviously prefer alternating cubic lattice. 
Furthermore, classical Molecular Dynamics which repro- 
duces their static correlation functions observed on the 
lattice (see discussion and references in [3) complement 
the Coulomb potential by the repulsive core, to ensure 
classical stability of the lattice. The potential is thus 
written as 



9rn 

47r 



(^)Q^+QJ 1 

h — 1 



(27) 



Se — m J ds — mj3 a/I -I- o? j fP- — Sc 



(25) 



where charges Qi, Qj = ±1 and additional dimensionless 
parameter b is the "core power" . It is usually selected 
to be large, to produce small corrections except close to 
the origin: we use the conventional b=9. The forces are 
balanced at distance 1, which defines the normalization. 

Near Tc the condensate is still small fraction of all 
particles, in spite of divergent (or highly peaked) spe- 
cific heat: thus we may think that only particles of one 
kind (e.g. only positive ones) are moving along the line. 
This prompts us to take one diagonal of the cubic lat- 
tice (rather than a line parallel to axes) to be jumping 
the distance v^, while keeping all other positive and all 
negative ones stationary. The potential corresponding to 



such forces was calculated and plotted in Fig. 14 (points) 
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FIG. 14: The effective potential for a diagonal moving line of 
charges (points) is compared to a fitted sinusoidal potential 
(curve). 



As usual, it is well represented by the sinusoidal potential 



(28) 



and therefor one can use the results obtained for ''He in 
Section |lT] (note that because of the diagonal supercur- 
rent we have substituted d \/2d into it). In particular, 
we will use caloron amplitude ( 20 1 . This formula was de- 
rived starting from a non-rclativistic Lagrangian, while 
we have shown that in this regime monopoles move rela- 
tivistically. In order to give an estimate of the relativistic 
correction we can compare the monopole mass predicted 
using Eq. [25] with the mass obtained from the analogous 
non-relativistic formula. 

Fixing Tc — 0.27 GeV and the jumping distance \f2d = 
we can extract the monopole mass as a function 
of the Feynman action. The result is presented in Fig. 15 
where we can see that the relativistic correction is very 
small, in the order of 5%, therefore we can hope that, 
also introducing the potential, using the non-relativistic 
formula the error in the prediction of the mass remains 
small. 

In order to include the potential we need an estimate 
of the coupling constant g^a- This estimate has been 
provided studying lattice correlation function in jl2l 113] 
to be such that plasma parameter 



r 



5; 



-^-KdTc 



(1-2) 



(29) 



In Fig. [16] you can see the dependence of the monopole 
mass from the Feynman parameter for three different 
choice of F: F = 0.5,1.0,2.0. Remember that for F ~ 1 
the plasma is a liquid (only for F << 1 we have a gas and 
for F ^ 100 a solid), therefore inside the range consid- 
ered the state of the system does not change, neverthe- 
less we can see that such a small variation in the plasma 
parameter is sufficient to produce a large variation in 
the monopole mass going from 50 MeV for F = 2 until 
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FIG. 15: Relativistic correction to the prediction of the 
monopole mass at T = Tc without the inclusion of effects 
connected with inter-particle coulomb potential. 



180 MeV for F — 0.5. On the other hand the dependence 
from the choice of the Feynman parameter is very small. 

In conclusion of this section we could say that in order 
to study the temperature dependence of the monopole 
mass it seems to be crucial to fix accurately the strength 
of the inter-particle potential. 
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FIG. 16: Monopole mass prediction at T = Tc using the non- 
relativistic formula for the Feynman parameter (Eq. 20 1. The 
three lines correspond to different strength of the coulomb 
potential. 



C. Discussion 

Fiquid "'He at the lambda-point shows characteristic 
infinitely high peak in its specific heat, famously at- 
tributed by Feynman to divergence in the sum over the 
contribution of the k-polygons. 

The pure gauge theory deconfinement is second or- 
der for SU(2) group, it is the first order transition for 
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SU{Nc > 2) and so far is designated to be a crossover one 
in QCD, with dynamical quarks. This means that added 
quarks somehow tame the peak in the specific heat, to 
still large but finite hight. (It may still be a finite- volume 
effect.) It is perhaps fair to say that at the moment no- 
body has a clue why is it so, at least we are not aware of 
any definite ideas. 

The real-world QCD has fundamental quarks, and for 
those another transition - chiral symmetry restoration - 
coincides or is very close to deconfinement . However for 
different quarks, with adjoint (roughly twice larger) color 
charge, those are separated and deconfinement is the first 
order transition. 

In this paper we assumed that deconfinement is the 
BEC of monopoles, and it is crucially important to test 
on the lattice whether it is indeed true or not. This 
means in general to look for monopole condensate at T < 
Tc- We would suggest to follow Feynman and testing his 
predictions: (i) at Tc the probability of polygons of any 
k gets comparable; and (ii) the extra action for "jumps" 
reaches the same critical value as for all other EEC's. 



We know from simulations for ^He that those were indeed 
true in this case: why should monopoles be any different? 
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Note added: When this paper was completed, we 
learned about one more interesting application of 
strongly coupled charged plasma possibly undergoing 
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